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5_^ , Abstract. Many phenomena in biology involve both reactions and chemotaxis. These pro- 

C^ ' cesses can clearly influence each other, and chemotaxis can play an important role in sustaining 

and speeding up the reaction. In continuation of our work |14) . we consider a model with a 
single density function involving diffusion, advection, chemotaxis, and absorbing reaction. The 
model is motivated, in particular, by the studies of coral broadcast spawning, where experi- 
mental observations of the efficiency of fertilization rates significantly exceed the data obtained 
from numerical models that do not take chemotaxis (attraction of sperm gametes by a chemical 
Ph ■ secreted by egg gametes) into account. We consider the case of the weakly coupled quadratic 

■^^ I reaction term, which is the most natural from the biological point of view and was left open in 

(-H ■ [M] . The result is that similarly to [M] , the chemotaxis plays a crucial role in ensuring efficiency 

of reaction. However, mathematically, the picture is quite different in the quadratic reaction 
case and is more subtle. The reaction is now complete even in the absence of chemotaxis, but 
the timescales are very different. Without chemotaxis, the reaction is very slow, especially 
for the weak reaction coupling. With chemotaxis, the timescale and efficiency of reaction are 
independent of the coupling parameter. 
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1. Introduction 



^ ■ Our goal in this paper is to study the effect chemotactic attraction may have on reproduction 

^ . processes in biology. A particular motivation for this study comes from the phenomenon of 
broadcast spawning. Broadcast spawning is a fertilization strategy used by various benthic 
invertebrates (sea urchins, anemones, corals) whereby males and females release sperm and egg 

k>l ' gametes into the surrounding flow. The gametes are positively buoyant, and rise to the surface 
^ . of the ocean. The sperm and egg are initially separated by the ambient water, and effective 

- - -' mixing is necessary for successful fertilization. The fertilized gametes form larva, which is 
negatively buoyant and tries to attach to the bottom of the ocean floor to start a new colony. 
For the coral spawning problem, field measurements of the fertilization rates are rarely below 
5%, and are often as high as 90% [3 dSl [221 ED] • On the other hand, numerical simulations based 
on the turbulent eddy diffusivity [B] predict fertilization rates of less than 1% due to the strong 
dilution of gametes. The turbulent eddy diffusivity approach involves two scalars that react and 
diffuse with the effective diffusivity taking the presence of the flow into account, instantaneous 
details of the advective transport not captured by eddy diffusivity approach. It is well known. 
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however, that the geometric structure of the fluid flow lost in the turbulent diffusivity approach 
can be important for improving the reaction rate (in the physical and engineering literature see 
[201 IMl EH]; in the mathematical literature see [I6l [3l HSl [9l |13] for further references). Recent 
work of Crimaldi, Hartford, Cadwell and Weiss [U [5] employed a more sophisticated model, 
taking into account the instantaneous details of the advective transport not captured by the 
eddy diffusivity approach. These papers showed that vortex stirring can generally enhance the 
reaction rate, perhaps accounting for some of the discrepancy between the numerical simulations 
and experiment. 

However, there is also experimental evidence that chemotaxis plays a role in coral fertilization: 
eggs release a chemical that attracts sperm PQ [21 [13 [IB] • Mathematically, chemotaxis has been 
extensively studied in the context of modeling mold and bacterial colonies. Since the original 
work of Patlak [2Tj and Keller-Segel [TDl E] where the first PDE model of chemotaxis was 
introduced, there has been an enormous amount of effort devoted to the possible blow up and 
regularity of solutions, as well as the asymptotic behavior and other properties (see [23] for 
further references). However, we are not aware of any rigorous or even computational work on 
the effects of chemotaxis for improved efficiency of biological reactions. 

In this paper, we continue the study of a simple single partial differential equation modeling 
the fertilization process that we initiated in [1^. The equation is given by 



9tp + u ■ Vp = Ap + xV ■ (pV[(A)-V]) - ep^ p(x, 0) = po(a;), x E R". (1.1) 

Here, in the simplest approximation, we consider just one density, p{x,t) > 0, corresponding 
to the assumption that the densities of sperm and egg gametes are identical. The vector field u 
in (II. ip models the ambient ocean fiow, is divergence free, regular and prescribed, independent 
of p. The second term on the right is the standard chemotactic term, in the same form as it 
appears in the (simplified) Keller-Segel equation (see [23]). This term describes the tendency 
of p(x, t) to move along the gradient of the chemical whose distribution is equal to — A~^p. 
This is an approximation to the full Keller-Segel system based on the assumption of chemical 
diffusion being much faster than diffusion of gamete densities. The term (— ep'') models the 
reaction (fertilization). The parameter e regulates the strength of the fertilization process. The 
value of e is small due to the fact that an egg gets fertilized only if a sperm attaches to a certain 
limited area on its surface (estimated to be about 1% of the total egg surface in, for example, 
sea urchins eggs [26]). We do not account for the product of the reaction - fertilized eggs - 
which drop out of the process. We are interested in the behavior of 

"^o(^) = / p{x,t)dx, 

which is the total fraction of the unfertilized eggs by time t. It is easy to see that mQ{t) is 
monotone decreasing. High efficiency fertilization corresponds to mQ{t) becoming small with 
time, as almost all egg gametes are fertilized. 

The case q > 2 in (II. ip has been considered in [13]. In this case, we proved that if chemotaxis 
is not present, there exists a constant fio{^,<l,d, po) such that mo(t) > po for all t. Moreover, if 
we fix everything except let e — )■ 0, po converges to mo(0). Hence, in the absence of chemotaxis 
weak reaction coupling leads to a very weak reaction completion rate. On the other hand, we 
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showed that ii d = 2 and chemotaxis is present, then mo(t) — )■ z/(x,g, po?^) as t — t- oo. Here, 
ly tends to zero if we increase x keeping everything else fixed. Moreover, u and the time scale 
of the convergence do not depend on e, showing drastic difference in behavior compared to the 
pure reaction-diffusion-advection case. 

The case d = q = 2 in (11. ip is the most natural one - it corresponds to the reaction propor- 
tional to the product of egg and sperm densities. It turns out that the situation in this case is 
more subtle. Without chemotaxis, it is not necessarily true that mQ{t) is bounded below from 
zero for all times by a fixed positive constant. Instead, mQ{t) — )■ as t — )■ 0, but very slowly. 
To formulate our results, let us define F{z) = z / log{{ez)^^) . 

Theorem 1.1. [Upper bound in the absence of chemotaxis] Suppose that p{x,t) satisfies 

dtp + u-Vp = Ap-ep\ p{x,0) = po{x) > 0, (1.2) 

whereu{x,t) G C°°(R^x [0, oo)) is divergence free, and its stream function H is globally bounded. 
Let the initial data po be such that 



/ po(a;)e"'^' dx < oo. 



for < a < ao- Then ||p(-,t)||ii — )■ as t — )■ oo. Moreover, there exists to > that depends 
only on a^ so that for all times t > t^ we have 

f(l|p(-.t)|k.)< ,^ „, fii'ii'r M'it'\i u,fV (1-3) 

l + eC(po)F(||p(-,to)||Li)log(t/to) 

Remarks. 1. The origin of the condition on stream function will be clear from the proof. An 
additional constant drift can be accommodated without any problem. Hence, the theorem 
holds, for example, for every smooth periodic fiow u. 

2. In particular, the bound (II. 3p implies that ||p(-, t)||2,i < C log log t/ log t for large t. 

3. The condition on the decay of po can be relaxed at the expense of making the bound in (II. Sp 
slightly weaker. For example, it will be clear from the proof that for po G 5 (Schwartz class), 
we have that for t > 1 and every a > 0, 

M-.OIU.< ,f';-"°L . (1.4) 

(1 + elogt) 
The bound (I1.3P is pretty close to being sharp, as the next theorem shows. 

Theorem 1.2. [A lower bound in the absence of chemotaxis] Suppose that p{x,t) satisfies (II. 2p 
where u{x,t) G C°°(R^ x [0, oo)) is divergence free. Suppose that the initial data po > G 5. 
Then there exists < to < oo depending on po and e such that for t >tQ we have 



||p(-,t)||Li > i + c'e||p(-,to)lui"log(tAo)exp(e||p(-,to)||LO' ^^'^^ 

Moreover, there exists Cq > so that z/e||po||Li < Cq, then 

\p{-,t) Li > — — - 1.6 

l + C;e||po||Lilog(t/to) 
fort > to = C*||po|lii/||po|li2, where C is a universal constant. 
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Hence the estimate of ~ (1 + elogt) ^ decay rate of L^ norm of p for large t is almost sharp. 
In the case when chemotaxis is present, we prove essentially the same result as in 



Theorem 1.3. [Estimates with chemotaxis] Let d = 2, and suppose that u G C°°(]R'' x [0, oo)) 
is divergence free. Assume that q = 2 and p{x,t) solves (11. II) with Po>OeS. Then 
a. If u = 0, then for every t > 0, we have 



IK.-)IU.<^(i + /i + ^). (1.7) 

b. If u ^ 0, then hm^^oo ||p(")^)||li < C{u,m2)x'^^^- Moreover, for < r < x^/^ we have 

\\p{;T)\\L.<Ciu,m,){xr)-'/\ (1.8) 



As in [13], the key feature of Theorem 11.31 is that all bounds are independent of e. Hence, 
chemotaxis provides a significant improvement over f ll.6p for reasonable time scales. The key 
ingredient in the proof of of this theorem in [13] was global regularity of solutions to fll.ip with 
d = 2, g > 2. It is well known that solutions of Keller-Segel equation can form singularities in 
finite time, but in our situation the reaction term has regularizing effect. However the proof in 
[13] relied on global in time control of L°° norm to prove global regularity, and this argument 
needed the condition g > 2 for reaction to offset the chemotactic nonlinearity. Here, we will 
prove global regularity for the case d = q = 2 using a different approach. Notice that the 
regularity properties of solutions to Keller-Segel type equations with reaction terms have been 
studied earlier in [251 I2SI EZ] , but in a different setting. 

2. The reaction-diffusion-advection case 

In this section, we prove Theorems 11.11 and 11.21 Let us begin with Theorem 11.11 In all 
arguments below C will denote universal constants that may change from line to line in the 
proof. Consider the advection-diffusion equation in M^ 

9t0 + M- V0 = A0, 0(x,O) = 0o(a;). (2.1) 

Suppose u is smooth and u = V^H, so that if is a stream function corresponding to u. Let us 
denote p{x,y,t) the fundamental solution of (12. ip : 



0(x,t) = / p{x,y,t)(j)o{y)dy. 

Recall the following well known bound on the fundamental solution of (12. ip . 

Theorem 2.1 (Norris [H]). Let(f){x,t) satisfy (12. ip . and assume that the stream function H{x) 
is bounded. Then 

p{x, y, t) < Ct'^ exp(-|x - y\^/Ct), (2.2) 

where C depends only on ||if||/,oo. 

The reference [19] contains more general results, but we stated the minimum that we need 
here. We have 
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Corollary 2.2. Suppose (po^x) > zs such that 

6°'^''^^ (j)o{x) dx < oo, 
for some < a < Uq. Then under conditions of Theorem \2.1\ we have 

f e"l^l0(x,t)dx<Ce^"'* / e"l^l0o(a;)rfa;. (2.3) 

Proof. A direct computation shows that 

e"^i0(a;,t) dx < Ce^"'* [ My)e"^' dy. 

JlR2 

Combining this and similar computations with —xi and ±X2 in the exponent leads to (12.3p . D 
Now suppose that p{x, t) solves 

^ip + M- Vp = Ap-ep^ (2.4) 

and 

/ po(x)e"''^' dx < oo, 

for every a < ao- 

Proposition 2.3. There exists to(ao) such that if t > to, then 



tlog(||p(-,t)||^J 

Proof. Consider 4>{x,t) solving (12.11) with initial data po{x) > 0. Then by comparison principle, 
p{x,t) < (j){x,t), and therefore by Corollary 12.21 we have 



^okl, 



p(x, t) dx < Ce^°'* / po(z/)e"l^l dy (2.6) 



for every « < ao? ^ ^ 0. We would like to find R{t) such that 



p{x,t)dx>- I p{x,t)dx, (2.7) 



where -BR(t) is a ball of radius R(t). Observe that by ( 12. 6p . 

p(x, t) rfa; < Ce^°'*-°^W / Po(y)e°l^l rfy. 






Thus choosing i?(t) so that 

i?(t)> Cat + a-Mog(||p(-,t) 11^,1) + C(po,«) 
will ensure (12.71) . Choosing 



a = t-'/W^ogiM;t)\\-') 
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optimizes the estimate provided that t > to(ao) so that a < oq. Then, for 

R{t)=Cipo)^Jtlog\\pi;t)\\ll 

we have 



Bnn) \Jb^,,. / \BBit)\ tlog{\\p{-,t)\\J] 



D 



Now we are ready to prove Theorem ll.li 
Proof of Theorem \l.l\ From Proposition I2.3[ we know that for all t >tQ, 

om-Ml^ = 4p{-Mh < -eC'(po)rip(-,t)iii,/iog(||p(-,t)||^,^). (2.8) 

Let us denote F{z) = z / (\og{ez)~^) . Then (12.81) implies 

f(iiM-.i)iu.)< ,^ „, ^'!l';i'V°'lii'\i ,,, ■ (2-9) 

l + eC(po)F(||p(-,to)i|Li)logV^o 

D 

Next, we prove Theorem II. 2[ 

Proof of Theorem M.B. Suppose p{x,t) satisfies (12.41) with po > E S, and (p{x,t) solves (12. ip 
with the same initial data. Then 

W(-.*)lli. = -llv^lli. <-"*-- "* 



lii llPollii 

by Nash inequality in dimension two and conservation of the integral of p. It follows that 

||0(-,t)||i.<min(||po||i.,Cl|po||i.ri). 

By duality (and incompressibility of u), we also have 

||0(-,i)||L- < min (||po||L°°,C'||po||L2i"^^^) • 
Combining these estimates via time split at t/2 yields 

||0(-,^)||l-' < min (||po||Loo,C||po||Lii"^) • 
By comparison principle, the same bound applies to p{x,t). Then we can estimate 

M-Mh < ||p(-,t)iu^||p(-,t)|Ui < c\\p{;t/2)U4p{-ML^r'- (2.10) 

A strictly weaker bound 

\\pi-Ml^<C\\po\\L4pi;t)\\L^t-' 
is also true. It implies that 

9,||p(-,t)|Ui = -e||p(-,t)||i. > -6C||po||Li||p(-,i)||Lit-\ 

and therefore 

||p(-,t/2)|Ui<2||p(-,t)|U.e^M.^ 
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Substituting this into (I2.10p . we get 

||p(-,t)||i. <Ce^"«'II^MIp(-,t)||i.t-^ 
Hence, we have 

OM-Ml^ > -6min(||po||i.,Ce^"^"ll-Mlp(-,i)lliii-')- 
Define to by 

Then, for t > to, we have 

leading to 

\\p{;to)\\L^ 



llp(-,t)IUi> 



1 + Ce||p(-, to) IIli log(t/to) exp(e||p(-, to) IIlO ■ 
If ellpollii < Co, with a sufficiently small Cq, then 

to < C||po||ii/||po||i2, 

and 

||p(-,to)||Li> ||po||Li(l-eC||po||Li). 
Hence, dLE]) follows. D 

3. The chemotactic case 

The proof of Theorem 11.31 is identical to the case q > 2 considered in [H] , provided that we 
have global regularity of solutions. It will be convenient for us to work in a setting similar to 
[Hj. Recall that we are considering the equation 

dtp + u-Vp = Ap + xV- (pV[(A)-V]) - ep', p{x, 0) = po(x), X G R\ (3.1) 

Define 

ii/iim„= / (ip(x)i + ivp(x)i)(i + |xr)dx, 

and the Banach space Kg^n is defined by the norm H/Ha'^,, = ||/||Af„ + ||/||ii'»- 

Theorem 3.1. Assume that n > and s > d/2 + 1 are integers and po > G Ks^n- Suppose 
that u G C°°(]R'^ X [0, oo)) is divergence free. Then there exists a unique solution p{x,t) of the 
equation (^ m C{K,^n, [0, oo)) n C°°(R"' x (0, oo)). 

The local existence is proved in the same manner as in [T3]. Also, it was shown in [13] that 
to prove global existence, it suffices to prove an a-priori global bound on Sobolev norms of the 
solution. 



8 



ALEXANDER KISELEV AND LENYA RYZHIK 



Proof. The first observation is that, integrating in space and time, we obtain a bound 

e / / p{x,tfdx < IIpoIIli, 

Jo JlR2 

for any smooth decaying solution of (13. ip on [0,T] x M^. This will be our main control which 
turns out to be sufficient for global regularity. Multiply equation (13. ip by (— A)*p and integrate. 
Let us denote by if* the homogenous Sobolev norm. We get 



^dMls<X 



[ (Vp)-(V(-A)-V)(-A)Vrfx +x f p'i-AYpdx +C{u 

JR2 Jm.2 

Consider the expression 

/ p\-Aypdx. 

Integrating by parts, we can represent this integral as a sum of terms of the form 

D^pD'-^pD'pdx, 
where / = 0, . . . , s and D denotes any partial derivative. By Holder inequality, we have 



2 II ||2 

H" ~ llPlliJs+i- 



D'pD'-'pD'pdx 



< IIpII^.IID'pIIl-.P P 



m-. 



with any pi and qi satisfying Pi^ + Qi ^ = 1/2. Observe that for d = 2, every integer < m < s, 
and 2 < p < oo, we have the Gagliardo-Nirenberg inequality 



in^/ill ^ /^ll /ill -i-^"!! /ill<^ 



m- - + 1 
p 



Then 

and therefore 

Next, consider 



s + 1 
\\D'p\\lpi\\D'-'p\\l^i <C||p|U2||p||^.+„ 

/ p'^{-Aypdx < C||p||i2||p||^,||p||^,+i. 

7r2 

/ (Vp)-(V(-A)-V)(-A)Vrfx. 

JR2 



Integrating by parts, we get terms that can be estimated similarly to the previous case, using 
the fact that the double Riesz transform dij{—A)~^ is bounded on L^, 1 < p < oo. The only 
exceptional terms that appear which have different structure are 

{d,, . . . di^p) ■ (V(-A)-V)5., . . . d.^pdx 



but these can be reduced to 



/ {di^...di^pYpdx 

Jr2 
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by another integration by parts, and estimated as before. Altogether, we get 

7;dt\\p\\]js < CCl|p||L2||p|lij«llp|lH«+i +Ciu)\\p\\jjs - ||p||^,+i < 



2 



ice\\p\\i^+c{ump\\%^+c{u)\\p\\i.-^ 



2 



Notice that 



B + d/2 



I nil <r llnll ''+l+''/2 llnlM+l+''''2 

\P\\h^ -^ WPWfjs+i WPWl'^ ' 



1 2 I ^ 

-^ - " "2 ^/'r<f2||^||2 I r^/„,\MI^I|2 I r</'„,MI^I|2 „ll ^11 s+i+d/2 



and so due to non-increase of L^ norm 

^^tiipiii. < (cewpwi^ + ciumpwis + c{u)\\p\\i. - c„,„^. 

From this differential inequality and integrability of ||p||^2 in time, a global bound for HpH^fs 
follows. n 
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